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Compressible isothermal magnetohydrodynamic turbulence is analyzed under the assumption of 
statistical homogeneity and in the asymptotic limit of large kinetic and magnetic Reynolds numbers. 
Following Kolmogorov we derive an exact relation for some two-point correlation functions which 
generalizes the expression recently found for hydrodynamics. We show that the magnetic field brings 
new source and flux terms into the dynamics which may act on the inertial range similarly as a source 
or a sink for the mean energy transfer rate. The introduction of a uniform magnetic field simplifies 
significantly the exact relation for which a simple phenomenology may be given. A prediction for 
axisymmetric energy spectra is eventually proposed. 
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I. INTRODUCTION 

Hydrodynamic turbulence, despite its ubiquitous na- 
ture, is extremely complex to be studied analytically. 
The degree of complexity gets considerably enhanced 
when the system consists of a magnetohydrodynamic 
(MHD) fluid. Yet several theoretical (analytical) works 
have been carried out in the framework of both the above- 
said cases when one assumes the incompressibility (den- 
sity of the fluid is constant in space and time and the 
velocity field is a solenoidal vector) of the fluid. It was 
Kolmogorov [1] who for the first time derived in 1941 an 
exact relation for the third-order moment of the veloc- 
ity structure functions in incompressible hydrodynamic 
turbulence - the famous 4/5 law [see also 2]. Later two 
important exact relations were derived respectively for 
the conduction of a passive scalar in a turbulent system 
[3] and incompressible MHD turbulence [4]. Note that 
exact relations can also be derived for quasi-geostrophic 
flows [5] or dispersive MHD [6, 7], whereas attempts are 
currently made to model axisymmetric turbulence [8]. 

In a recent Letter [9] we have derived, for the first time, 
an exact relation for compressible hydrodynamic turbu- 
lence with an isothermal analytic closure P ~ c^p, where 
P, p and Cs are respectively the fluid pressure, density and 
sound speed (for our case it is supposed to be constant) 
defined at every point of the flow field. In the present ar- 
ticle we generalize our previous work to the MHD case. 
Of course, the taking into consideration of compressibil- 
ity of an MHD fluid leads us to get closer to the reality 
taking place for example in the solar wind because (i) it is 
a plasma (so contains charged species), (ii) Ulysses data 
analysis has shown clear evidences of the effect of com- 
pressibility [10] and (iii) Voyager data analysis at several 
Astronomical Units (AU) reveals the presence of both 
fluctuations and jumps with spectra steeper than Kol- 
mogorov [11]. Moreover turbulence inside the interstellar 
clouds, being highly compressible (supersonic), also de- 
mands a theoretical background in order to be properly 



understood [12-14]. A thorough theoretical work in com- 
pressible MHD turbulence is essential (suggested in [15]) 
above all because it will help to analyse and understand 
the results of Direct Numerical Simulations (DNS) done 
for compressible turbulence [16-20]. 

A number of numerical works have been carried out for 
the last two decades in the case of compressible hydro- 
dynamic turbulence. Being initiated mainly around the 
years 1990 [21, 22], the numerical simulations get con- 
siderably sophisticated in the last ten years. Recently, 
Kritsuk et al. [18] have studied three-dimensional (3D) 
isothermal compressible hydrodynamic turbulence using 
Piecewise Parabohc Method (PPM) for upto 2048^ grid 
points where the r.m.s. Mach number is around 6. Their 
study consists of the Helmholtz decomposition of veloc- 
ity and they obtain a velocity power spectrum around 
-2 (-1.95, -1.92, -2.02 respectively for total velocity, 
solenoidal and compressible velocity components) using 
1024'^ grid points whereas they have a —5/3 spectrum for 
the density-weighted fluid velocity p^^^v. Another impor- 
tant numerical work has been accomplished by Federrath 
et al. [19] for an identical (to the previous one) system of 
turbulence using upto 1024^ grid points by applying com- 
pressible or solenoidal forcing. For two different types of 
forcing, they have obtained remarkably different power 
spectra (which questions the universality) corresponding 
to various quantities like the total velocity, its longitu- 
dinal and transverse components, Inp, p^/^w or p^^^v. 
Total velocity gives a —1.86 and a —1.94 power spectra 
in the supersonic range for solenoidal and compressible 
forcings respectively which are close to the Burgers spec- 
trum (in k~^). The compensated spectra in the subsonic 
region cannot be affirmed to date to be free from bot- 
tleneck effect and so it is hard to conclude anything for 
that part. However, a —5/3 spectra is obtained for p^^^v 
even in supersonic range using solenoidal forcing whereas 
a compressible forcing gives a —2.1 spectrum. Here too 
a plausible presence of bottleneck effect prevents us from 
inferring for subsonic region. Recently, Aluie et al. [23] 
have provided empirical evidences of kinetic energy cas- 
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cade (which is not an invariant in compressible turbu- 
lence) using subsonic and transonic simulations. Such 
cascading takes place far in the incrtial zone due to very 
fast decay of the pressure-dilatation term which domi- 
nates at large scales. They do not give, however, any 
intuition for supersonic turbulence for which a priori the 
kinetic energy cascade may not be the dominant process. 
In the meantime, a theoretical work based on coarse- 
graining [24] shows that for the kinetic energy a range 
of scales independent of the dissipation and the forcing 
may exist in compressible turbulence. 

3D DNS of compressible MHD turbulence is more dif- 
ficult to perform since it requires the computation of a 
heavier system in which several new regimes are possible. 
In this framework, for example Cho and Lazarian [25] im- 
plemented a third-order accurate hybrid and essentially 
non oscillatory numerical scheme to solve the isothermal, 
ideal MHD equations in a periodic box and for obtaining 
scaling relations they drove turbulence with a solenoidal 
forcing in Fourier space using 256^ grid points. In their 
study, they executed a mode separation method in order 
to study the individual scaling behaviour of each of the 
present linear modes in a compressible MHD system i.e. 
Alfven mode, slow and fast magnetosonic (MS) modes. 
They reported a —5/3 velocity fluctuation spectrum for 
each of the Alfven and slow MS modes whereas a —3/2 
spectrum for fast MS mode. The nature of spectra ob- 
tained is found to depend very weakly on the value of the 
plasma f3 parameter (the ratio of kinetic pressure to mag- 
netic pressure). An approximately {£ is the length 
scale) scaling is also obtained for the second-order struc- 
ture function. In [26] a Kolmogorov spectrum (—5/3) was 
reported for both the velocity (v) and density-weighted 
velocity (p^^^v) fluctuations for subsonic MHD turbu- 
lence whereas only the second quantity follows a —5/3 
spectrum for supersonic turbulence. 

A relevant work to our present paper is that of Kritsuk 
et al. [27] where the authors studied numerically super- 
sonic turbulence in magnetized magnetic clouds. The 
work is realised with grid resolutions from 256'^ upto 
1024'^ using PPM, an isothermal closure and a solenoidal 
forcing to drive the fluid. Interestingly, they have ob- 
served the scale invariance of the total (kinetic -I- com- 
pressible) energy transfer rate in the inertial zone; we 
shall use this concept later in this article for construct- 
ing the different spectra. The most interesting and im- 
portant aspect of their article is to obtain numerically 
the linear scaling in the 4/3-law of incompressible MHD 
turbulence to the weakly magnetized supersonic turbu- 
lence - the former Elsasser variables (z* ) being replaced 



tion HI; Section IV is devoted to a discussion about the 
main result where in particular we look at the effect of 
a uniform magnetic field; finally a conclusion is drawn in 
Section V. 



II. COMPRESSIBLE MHD EQUATIONS 

For constituting an exact relation of compressible 
MHD turbulence of an isothermal plasma, we rewrite the 
basic equations (except the continuity equation which is 
kept unchanged) of ideal MHD in an unusual form [28] 
as follows: 



dtp + V-{pv) = 0, 



(1) 



1 / 1 , 

VA • V VA V P + -pVA' 

P \ 2 



-va(V • va) + dK 

VA • V V — 

-2p(V • va) , (4) 



fo , (2) 
KV-v)+dM, (3) 



dtv + V • V V = 

i9tVA + V • V VA = 
VA • V p = 

with p the plasma density, v the fluid velocity field, 
Va = b/^^oP the Alfven velocity and P the pressure. 
The terms dK , dM and fo represent respectively the con- 
tribution of the kinetic viscosity, magnetic resistivity and 
the external forcing. It can be identified that equations 
(2), (3) and (4) represent the force equation, the Fara- 
day's equation (with idealised Ohm's law injected) and 
the zero divergence of magnetic field respectively. For the 
sake of simplicity in the present paper we use the isother- 
mal closure {i.e. P = C^p). It is important to note that 
va as it is defined in our case does not correspond to the 
usual Alfven velocity which does not depend on the den- 
sity but just on the mean density. Moreover, the Alfven 
mode is defined only when there is a mean external mag- 
netic field Bo whereas our va is defined independently 
of the presence of Bo. 

As in the hydrodynamic case, our analysis will be car- 
ried out in the physical space in terms of two-point cor- 
relation functions and structure functions, where the un- 
primed c^uantities represent the properties at the point 
X and the primed quantities correspond to the point x' 
(with x' = x + r). Our system is supposed to be statis- 
tically homogeneous and to undergo a completely devel- 
opped turbulence. The analysis is general and contains 
no hypothesis of isotropy. 

In order to establish an exact relation, we shall utilise 
the conservation of the total energy whose density is 



by the density- weighted Elsasser variables (p^^'^z^). The given as: 
present article will give, by the help of an exact relation 
for compressible isothermal MHD turbulence, a theoret- 
ical support for this generalization. 

The rest of the paper is organized as follows. Section 
II introduces the compressible MHD equations and the 
correlation functions; the derivation of the exact relation 
which is the main result of this paper is given in Sec- 



£^(x) = -(v • V + VA ■ Va) 



pe, 



(5) 



where e accounts for the compressible energy term due 
to the isothermal closure and is expressed as ln(p/po), 
where po is the initial equilibrium density. We however 
introduce another quantity which wc call compressible 



3 



cross-hclicity density and is written as: 

i7(x) = pv • VA . (6) 

Although this quantity is not conserved in compressible 
MHD, we need it (as we shall see later) to express the 
final flux term in a more indicative way (in connection 
with the incompressible exact relations [4] ) . We can write 
the expressions for E'{x') and -ff'(x') in an identical way 
as that of the above. 



III. DERIVATION OF AN EXACT RELATION 

Following the formalism of our previous paper [9] , we 
define the relevant two-point correlation functions asso- 
ciated to the total energy and the compressible cross- 
helicity density as: 

Re = ^(v-v' + VA-VA')+pe', (7) 
Rh = ^(v-va'+va-v'), (8) 

and similarly for the primed quantities. Using the above 
definitions and introducing compressible Elsasser vari- 

I 



9t(pv • v') = 



i9t(pvA • va') = 



dt{pv ■ Va') = 



dtip^A ■ v') = 



dtipe') = 



where di and fi{i = l to 4) represent the dissipative and 
forcing contributions respectively. Hence we have all the 



ables = V zL va, we obtain: 

{E + E')±{H + H') - {Re ± Rh + Re ± Rh) 

-SpSe = i(5(pz±) • Jz± , (9) 

where for any variable ip, S^p = Tp{'x.-\-r) — ipix) = ip' — ijj. 
Under ensemble average (which is equivalent to spatial 
average for a statistically homogeneous system) the above 
relations get reduced to: 

{E)±{H) - ^{RE±RH + R'E±R'H)-l{SpSe) 

= l(%z±).5z±). (10) 

Now we shall calculate: 

dt{RE ± Rh + R'e ± R'h) , (11) 

in order to see how the incompressible relations [4] can 
be generalized to the compressible pseudo-energies (e^ = 
pz^ ■ z*) which are no more conserved. By using the 
equations (l)-(4) and basic identities of vector calculus, 
we obtain after some calculations: 



I 

necessary elements for the dynamic equation of Re and 
Rh- Similarly we can write (by symmetry) the expres- 



- V • [p(v • v')v] + V • [pvA(v' • va)] - (V • va) [V • (pva)] - V • (Pv') - V • (-pvA^v') (12) 

- p(v' • va)(V • va) - V • [p(v.v')v'] + p(v • v')(V' • v') + V • [p(v • va')va'] - p(v • va')(V' • va') 

- 4V' • (P'v) - 4v'(i Aa"v) - p(v • va')(V' • va') + + a , 

- V • [p(va' • va)v] + V • [pva(v • va')] - (va' • v)[V • (/ova)] - ^^(va • va')(V • v) (13) 

- ^p(vA • va')(V' • v') - V • [p(vA • va')v'] + p(vA • va')(V' • v') + V • [p(vA • v')va'] 

- p(vA-v')(V'-VA') + d2, 

- V • [p^r{v ■ va')] + V • [pva(va • va')] - (va' • va)[V • (pva)] - V • (Pva') - V • (^PW^va') (14) 

- p(va' • va)(V • va) - V • [p(v • va')v'] + p(v • va')(V' • v') + V • [p(v • v')va'] 

- p(v • v')(V' • va') - ^p(v • va')(V' • v') + d3 + /2 , 

- V • [/9v(v' • va)] + V • [pva(v' • v)] - (v • v')[V • (/ova)] - ip(v' • va)(V • v) (15) 

- V • [p(vA • v')v'] + (V • v')(pvA • v') + V • [p(vA • va')va'] - p(vA • va')(V' • va') 

- 4^' • [^'^a] - • [Ip'v'I^a] - P(VA • va')(V' • va') + d4 + /3 , 
P P 2 

- V • (pe'V) - V • (pe'v) - V' • (Pv') + pe'(V' • v') . (16) 
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sions of R'^ and R'^ just by carefully interchanging the tistical averaging operator on each term, after judicious 
primed and the unprimed quantites. Now if we apply sta- re-arrangement of all the terms we finally obtain: 



dt{RE±RH + i?B±i?'ff)=V^ 



1 



(5(/9z=^)(5z=^ + 6pSe 



(5zT+(5(e + ^),5(pz±) 



i(^V'.(pzV) + ^V.(p'z'V 



+ ^(V-v) 
+ ((V'-v') 
'(V-va) 
+ ('(V'-vaO 



R^^±R'^-E'TH'T^{i^A-z'^)-^ + ^ 



i?i,±i?ff -i?Ti?T ^(va' -z^) - I + ^ 

— 5P' P' 

Rh±Re-R'hTR'e±E' + H' - (5p(vA • z'±) ± — ± ^ 



R'h±R'e-RhTRe±E + H - d>(vA' ■z^)±^±'-^ 



5P . 



M 



(17) 



D±+F± 



where (5X ee (X + X')/2 and Pm = {l/'2)pv\ is the 
magnetic pressure. D± and F± represent respectively 
the averaged dissipative and forcing terms. The above 
expression consists of two equations which represent the 
generalized form (to compressible MHD) of the coupled 
equations (4) and (6) of [4]. As it was already shown 
[28] , for compressible flows the cross-helicity is no longer 
conserved and we cannot (by equation (6)) let the term 



dt{RE ± Rh + R'e i R'h) vanish to the extent of a sta- 
tionary state where the average forcing term and the 
average resultant dissipative (kinetic -I- magnetic) terms 
cancel each other to ensure the conservation of total en- 
ergy. Hence for the compressible case, these two equa- 
tions individually do not lead to an exact relation com- 
prising of a universal quantity. However, if we add both 
we find: 



dt {Re + r'e 



-(5(pz ) ■ Sz + SpSe 



l^ipz^ 



SpSe 

2 



Sz- + 6{e + -f)S{pz- 



pz' 



(18) 



1 Mv' • (pz+e') + 4v • (p'z'+e) + %V ■ {pz- e') + %V • {p'z!-e\ 

o\C„ C„ C^. C„ 



+ ((V-v) 



i?R — E' 



5p, , . P'm-P' 
^(va'-va) + ^^ 



+ (V'-v') 



RE-E-—{yfj^-WA) + ^ 



((V • va) [Rh -R'h + H' - 5p{w' ■ va)] ) + ((V • va') [R'h-Rh + H- 5p{w ■ va')] ) + D + F, 



where D = {D+ +D-)/2 and F ^ {F+ + F_)/2. Now 
to establish the final relation, we assume that in the limit 
of infinite (kinetic and magnetic) Reynold's numbers, the 
system attains a stationary state. This assumption, when 
applied to the relation (18), shows that the left hand term 



of the above equation vanishes. Now if we concentrate 
far in the inertial zone where the dissipative terms are 
negligible with respect to the other terms, we are left 
with: 
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-2£ 



-(5(/9z ) ■ 5z + 5p5e 



-(5(pz^) • 5z^ + SpSe 



Sz- 



i • (pz+e') + • (p'z'+e) + • (pz-e') + • (p'z'-g) 



Sie+^)5{pz 



pz' 



(19) 



(V-v) 



i?'^ — E' 



I ^ P'm-P' 



+ ((V'.v') 



i?B - - y (VA • VA ) + ^ 



((V • va) [Rh - R'h + H' - 5piV ■ va)] ) + ((V • va') [R'h-Rh+H- <5p(v • va')] ) , 



where e denotes the mean total energy injection rate 
(which is equal to the mean total energy dissipation rate). 
In the above equation the flux terms arc deliberately 
written in terms of the compressible Elsasser variables z^ 
whereas the source terms are expressed in terms of the ve- 
locity and the compressible Alfven velocity. This writing 
is expected to be useful for (i) an attempt to generalize 
the Alfven effect (introduced by Kraichnan) to describe 
the phenomenology of compressible MHD turbulence and 
(ii) understanding separately the contribution of the ve- 
locity field and the Alfven term in the source term. It is 
interesting to note that in the second part of our flux, the 
/3 parameter {i.e. c'j./v\) appears which can help us un- 
derstand (or obtain) the different limits depending upon 
its value. The above equation, being considerably bulky, 
is laborious to handle. For quick references, we give be- 
low a schematic view of that equation: 



-2e 



$1 



1 



(V • v)5i + (V • va)52 



(20) 



IV. DISCUSSION 



of V and b (assuming the constant density to be nor- 
malised to unity) help us recognise the primitive form of 
the equations found in [4]. 



2. Compressible hydrodynamics 



To verify the hydrodynamic limit, we put va = and 
so z"'' = z 



-2e = 

+ /(V-v) 



v; we are left with: 



— (5(pv) • (5v + SpSe 



(5v + 6e5{p'v) 



P' 

R'e - P' - — 



(V'.v') 



Re — E 



P 
"2 



(22) 



It is nothing but the exact relation which was derived 
in [9] . The only difference between the relation obtained 
above and that in the published paper is that here the 
pressure terms are written as the source terms whereas 
previously those were considered to contribute in flux 
terms. 



A. Relevant limits 

Equation (19) is the main result of the paper. It is an 
exact relation for compressible isothermal MHD turbu- 
lence in terms of correlation functions. From this equa- 
tion we can immediately investigate particular limits. 

1. Incompressible MHD 

In this limit, all the source terms vanish. Moreover, the 
compressible energy term becomes zero (evident from its 
definition as for incompressible case p ~ po) and we get 
finally (normalising the density to unity): 

-2e = • {{6z-f6z+ + {5z+f5z-) . (21) 

Note that the term • {6{v\/2)5{z^ + z'^)) can be 
shown to be zero in incompressible case by expanding 
it and letting the enter inside the average operator 
in a judicious manner. A rewriting of the in terms 



3. High and low /3 plasmas 

Without problem we admit that in the limit where 
the beta parameter of the plasma tends to infinity (very 
large value) , i. e. the plasma becomes almost incompress- 
ible (although not entirely), the flux term 4>2//3 becomes 
negligible with respect to $i of equation (20). On the 
contrary for a very small beta value, where the plasma 
can be assumed to be cold and magnetised (kinetic pres- 
sure negligible with respect to magnetic pressure), the 
term <I>2//3 dominates over $i term and at that situation 
the effective flux term becomes (after rearrangement): 

-^(^V'.(pve') + iv.(A'e)) . (23) 

B. Presence of an external magnetic field 

Relation (19) comprises the total magnetic field at each 
point of the fiow field. This total field b at each point can 
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be supposed (a realistic case) to have a fluctuating part 
(vary in space and time) b superimposed on a constant 
external magnetic field Bq. In the following, we shall 
investigate the flux and the source terms under the said 
situation. 



1. Flux contribution 



The part of the total flux term which contains Bq can 
be expressed as (with fiQ = 1): 



2 



1 



S ( j 6{^) [Bo X {6v X Bo)] + S{^) [Bq x {5v x Sva)] -S^i^-^] [S{pv) ■ Bq] Bo (24) 



1 



5 (^-^^ [(5(pvA) x ((5v X Bo) + 6{pv) x (5va x Bq)] + B^S Q ) 6(p^) + 2 



Bo-b ^, , , , ^ Bo b' , , „ \ 

—;^pv ■ Vp' + ^p'v' ■ Vp + —^p'^' ■ Vp 



B...fl|) 



5(pv) 



r 



where va = h/^/p. Now we assume the external field Bo 
to be very strong so that Bq |b| (and also Bq ^ \S-v\). 
This situation is classical in astrophysics: for example, 
there are many evidences of magnetic loops in the tur- 
bulent solar corona which are interpreted as strong uni- 
form magnetic fields on which small magnetic and ve- 
locity fluctuations are present [29, 30]. Wc shaU just 
consider the terms weighted by Bq . After some straight- 
forward calculations, we obtain the resultant fiux term 
(the magnetic terms without Bq and with single power 
of Bq are neglected) which writes at main order: 



2 



-TP 



1 



(5(Vp) [Bq X (5v X Bo)] 



[S{pv) ■ Bo] Bo ) . (25) 



The above expression gives the modifying part of the fiux 
in the presence of a strong constant magnetic field applied 
externally. One can easily understand that the modifica- 
tion is purely due to compressibility. (Note that the pure 
kinetic terms (not shown in (25)) gives also a contribu- 
tion to the total flux.) It is also interesting to notice that 
in expression (25) the fluctuations are exclusively kinetic 
in nature (because of the absence of e.g. Hall type term 
in the basic equations) . One can easily verify that the Sv 
of the first term and the S{pv) of the second term of the 
above expression can be replaced by ^v_l and ^ (pv|| ) re- 
spectively where (5vx_LBo and 5v|| ||Bo. The pure kinetic 
terms can however be omitted by assuming the external 
magnetic contribution to be dominant with respect to the 
velocity and the density fluctuations and then expression 
(25) will represent the total flux contribution. 



2. Source contribution 

The source terms are also modified due to the effect of 
a strong external magnetic field. At main order (keeping 



only the terms in Bq), the terms of type ((V • v)^!) get 
reduced to: 



B; 



1 



2^^<5(V.v)5^— )<5(Vp)-5(V.v)), (26) 



and the source terms like ((V • va)'S'2) writes: 



(27) 



Bo.(V(- 



(Bo.v')>'.(-L)}-(Bo.v)^; 



(B 



-(Bo-V) 



We note that the latter expression implies only parallel 
components of the velocity. 



3. Reduced anisotropic law 

Further simplifications are possible if we assume that 
the velocity field vector at each point of the flow field is 
(at the main order) perpendicular to the external mag- 
netic field i.e. if v\\ — = 0, and therefore v = v^. 
In that case, (4'va) vanishes and so is the second term 
of i^Bz)] then the corresponding reduced exact relation 
can be simply written as: 



-4e = B^Wr^ 



1 



S[ — ] 6{^)5v^ ) (28) 




5)..V^v,)(,.^^j), 

where V_l implies derivatives transverse to Bo. Expres- 
sion (28) is the second main result of this paper: it is 
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a limiting case of the exact relation (19) for isothermal 
compressible MHD turbulence under the influence of a 
strong external uniform magnetic field (sub-alfvenic tur- 
bulence regime). We may simplify the previous equation 
by assuming axisymmetry; the exact relation (28) can be 
written symbolically as (by using cylindrical coordinates) 



-4e 



1 



(29) 



where the radial component of the energy flux vec- 
tor (first term in the right hand side of (28)) and S is 
a source/sink term (last term in the right hand side of 
(28)). If we define an effective mean total energy injec- 
tion rate as Eoff = e -I- 5/4, a simple interpretation of 
expression (29) can be proposed as we see in Fig. 1: 
whereas for a direct cascade the energy flux vectors are 
oriented towards the axis of the cylinder, dilatation and 
compression are additional effects which act respectively 
in the opposite or in the same direction as the flux vectors 
(since terms like, 1 + p' / p, are positive). 



4. Phenomenology and energy spectra 

In this section we shall make a prediction on compress- 
ible spectra using the derived theoretical relation (28): 
although simplified, it is considerably indicative for phc- 
nomcnological intuition. 

The existence of strong external magnetic field Bq ren- 
ders the energy transfer in the parallel (along Bq) di- 
rection negligible in comparison to the transverse one 
[31-36]. Then, our prediction will be made for trans- 
verse spectra {i.e. with a spectral dependence only in 
ki_). From the flux term of expression (28), we find di- 
mensionally (keeping the source terms aside): 



(30) 



The above expression corresponds to a compressible phe- 
nomenology in which the magnetic energy density, piv\^, 
is transferred through the scales at the transfer time, 
Ti ~ and the transfer rate e. We may define the 

weighted variables: 



Bi = pe ' VAi = pf 
which leads to: 



1/6 



Bq , 



Whence, the spectral prediction: 

-5/2 



E''ik^)^E^{k^)^ek] 



(31) 
(32) 

(33) 
(34) 



Taking into consideration the source terms will only 
modify e for giving an effective mean total energy in- 
jection rate Ecff (likewise in compressible hydrodynamics 



n n n 




Seff<S 



Seff>S 



FIG. 1: Dilatation (left) and compression (right) phases in 
space correlation for axisymmetric MHD turbulence. In a 
direct cascade scenario the flux vectors (dotted arrows) are 
oriented towards the axis of the cylinder. Dilatation and 
compression (solid arrows) are additional effects which act 
respectively in the opposite or in the same direction as the 
flux vectors. 



[9]). A compressional effect (V • v < 0) will increase 
the effective energy flux rate whereas a dilatational ef- 
fect (V • V > 0) will reduce SeS as it is evident from 
the source terms of expression (28). As stated before a 
simple picture emerges (see Fig. 1) with the additional 
hypothesis of axisymmetry (symmetry in the plane per- 
pendicular to Bq). The flux vectors (dotted arrows) are 
always oriented towards the axis of the cylinder since a 
direct cascade is expected. Dilatation and compression 
act additionally (solid arrows): in the first case, the ef- 
fect is similar to a decrease of the local mean total energy 
transfer rate whereas in the second case it is similar to 
an increase of the local mean total energy transfer rate. 
Then, the final prediction for compressible turbulence un- 
der a strong uniform magnetic field is: 



E' 



\k^_)^jE^{k^) ^ e,s{k^)k 



-5/2 



(35) 



with a priori a possible dependence of the effective mean 
total energy transfer rate on k±. It that case, a power 
law steeper than —5/2 may be observed at large scales 
when compressible MHD turbulence becomes supersonic. 



V. CONCLUSION 

The present paper is a generalization of a recent work 
made for compressible hydrodynamics [9] in which it was 
possible to give a simple interpretation of the different 
terms of the exact relation. The introduction of both 
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the Lorentz force and the induction equation renders the 
situation more difhcult to analyze but it is beheved that 
the main physical characteristic of compressible MHD 
emerges when a uniform magnetic field is included. In 
that case, we have derived rigorously a reduced expres- 
sion (28) which can be used easily either to analyze direct 
numerical simulations or in situ observations like those 
obtained in the solar wind. Since most of the natural 
plasmas are in interaction with a large scale magnetic 
field, it is likely that the reduced form is the most inter- 
esting one. 

Our work confirms the relevance of density-weighted 
Elsasser variables, p^^^z^, for compressible MHD turbu- 
lence. However, unlike the hydrodynamic case, here one 
should not replace by p^/'^z^ for recovering the in- 
compressible scaling laws. Instead the appropriate vari- 
able should be, p^/^(z+^z^ +z^^z+)^/'^, for obtaining a 
satisfactory energy scahng as suggested by relation (19). 
The said relation holds good for any type of compressible 
MHD turbulence (subsonic, transonic, supersonic, sub- 
Alfvenic or super Alfvenic). Hence it can directly be used 
to analyze direct numerical simulations in compressible 



MHD turbulence at different Mach numbers or for plas- 
mas with different /3 parameter in order to understand 
the underlying phenomenology. In order to approach 
the astrophysical turbulence more efficiently we have to 
take into account the gravitational force and also a poly- 
tropic closure although they would increase enormously 
the complexity of the analysis. 

It was mentioned in the past that constant flux so- 
lutions are in general not relevant for compressible tur- 
bulence because the coupling to a sonic field provides a 
supplementary sink which will modify the inertial index 
[37-39]. A modification of the solenoidal energy spec- 
trum (for hydrodynamics turbulence) was proposed by 
introducing a simple ad hoc model in which the main 
compressible parameter was the Mach number. Here, we 
may reach a similar conclusion with a rigorous analysis in 
terms of two-point correlation functions. We believe that 
our interpretation in terms of effective energy flux rate, 
source and flux terms are the key ingredients to further 
understanding of compressible turbulence. 
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